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INFLUENCE OF BOTTOM TOPOGRAPHY ON LONG WATER WAVES 

Florent Chazel 1 

Abstract. We focus here on the water waves problem for uneven bottoms in the long- wave regime, on 
an unbounded two or three-dimensional domain. In order to derive asymptotic models for this problem, 
we consider two different regimes of bottom topography, one for small variations in amplitude, and one 
for strong variations. Starting from the Zakharov formulation of this problem, we rigorously compute 
the asymptotic expansion of the involved Dirichlet-Neumann operator. Then, following the global 
strategy introduced by Bona, Colin and Lannes in [6] , new symetric asymptotic models are derived for 
each regime of bottom topography. Solutions of these systems are proved to give good approximations 
of solutions of the water waves problem. These results hold for solutions that evanesce at infinity as 
well as for spatially periodic ones. 

Resume. Nous nous interessons ici au probleme d'Euler surface libre pour des fonds non plats en 
regime d'ondes longues, sur un domaine non borne a deux ou trois dimensions. Ann de construire 
des modeles asymptotiques pour ce probleme, nous considerons deux regimes topographiques sur le 
fond du domaine, l'un pour de petites variations en amplitude, et l'autre pour de fortes variations. A 
partir de la formulation de Zakhzarov, nous contruisons rigoureusement le developpement asymptotique 
de l'operateur de Dirichlet-Neumann relatif au probleme. En suivant la strategie globale proposee par 
Bona, Colin et Lannes dans [6] , nous obtenons ensuite de nouveaux modeles asymptotiques symetriques 
pour chaque regime de variation topographique du fond. Nous prouvons alors que les solutions de ces 
systemes fournissent de bonnes approximations aux solutions des equations d'Euler surface libre. Ces 
resultats sont valables aussi bien pour des solutions evanescentes a l'infini que pour des solutions 
spatialement periodiques. 
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Introduction 

Generalities 

This paper deals with the water waves problem for uneven bottoms which consists in describing the motion 
of the free surface and the evolution of the velocity field of a layer of fluid, under the following assumptions : 
the fluid is ideal, incompressible, irrotationnal, and under the only influence of gravity. 

Earlier works have set a good theoretical background for this problem : its well-posedness has been discussed 
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among others by Nalimov ( [18], 1974), Yoshihara ( [29], 1982), Craig ( [9], 1985), Wu ( [27], 1997, [28], 1999) 
and Lannes ( [16], 2005). 

Nevertheless, the solutions of these equations are very difficult to describe, because of the complexity of these 
equations. At this point, a classical method is to choose an asymptotic regime, in which we look for approximate 
models and hence for approximate solutions. We consider in this paper the so-called long-wave regime, where 
the ratio of the typical amplitude of the waves over the mean depth and the ratio of the square of the mean 
depth over the square of the typical wave-length are both neglictible in front of 1 and of the same order. 
In 2002, Bona, Chen and Saut constructed in [5] a large class of systems for this regime and performed a formal 
study in the two-dimensional case. A significant step forward has been made in 2005 by Bona, Colin and 
Lannes in [6]; they rigorously justified the systems of Bona, Chen and Saut, and derived a new specific class of 
symmetric systems. Solutions of these systems are proved to tend to solutions of the water waves problem on a 
long time scale, as the amplitude becomes small and the wavelength large. Thanks to their symmetric structure, 
computing solutions of such systems is significantly easier than computing directly solutions of the water waves 
problem. Another significant work in this field is the one of Lannes and Saut ( [17], 2006) on weakly transverse 
Boussinesq systems. 

However, all these results only hold for flat bottoms. The case of uneven bottoms has been less investigated ; 
some of the significant references are Peregrine ( [23], 1967), Madsen et al. ( [19], 1991), Nwogu ( [22], 1993), 
and Chen ( [8], 2004). Peregrine was the first one to formulate the classical Boussinesq equations for waves in 
shallow water with variable depth on a three-dimensionnal domain. Following this work, Madsen et al. and 
Nwogu derived new Boussinesq-like systems for uneven bottoms with improved linear dispersion properties. 
Recently, Chen performed a formal study of the water waves problem for uneven bottoms with small variations 
in amplitude, in ID of surface, and derived a class of asymptotic models inspired by the work of Bona, Chen 
and Saut. To our knowledge, the only rigorously justified result on the uneven bottoms case is the work of 
Iguchi ( [12], 2004), who provided a rigorous approximation via a system of KdV-likc equations, in the case of 
a slowly varying bottom. 

The main idea of our paper is to reconsider the water waves problem for uneven bottoms in the angle shown 
by Bona, Colin and Lannes. Moreover, our goal is to consider two different types of bottoms : bottoms with 
small variations in amplitude, and bottoms with strong variations in amplitude. To this end, we introduce a 
new parameter to characterize the shape of the bottom. In the end, new asymptotic models are derived, studied 
and rigorously justified under the assumption that long time solutions to the water waves equations exist. 

Presentation and formulation of the problem 

In this paper, we work indifferently in two or three dimensions. Let us denote by X £ ~R d the transverse 
variable, d being equal to 1 or 2. In the two-dimensional case, d — 1 and X — x corresponds to the coordinate 
along the primary direction of propagation whilst in the three-dimensional case, d = 2 and X = (x, y) represents 
the horizontal variables. We restrict our study to the case where the free surface and the bottom can be described 
by the graph of two functions (t, X) — > r](t, X) and X — > b(X) defined respectively over the surface z = and 
the mean depth z = —ho both at the steady state, t corresponding to the time variable. The time-dependant 
domain fl t of the fluid is thus taken of the form : 



In order to avoid some special physical cases such as the presence of islands or beaches, we set a condition of 
minimal water depth : there exists a strictly positive constant h m i n such that 



Q t = {(X,z), X£R d , 



-h + b(X)<z<r](t,X)} . 



7](t, X) + h - b(X) > h min , (t, X) £ K X K 2 . 



(0.1) 



For the sake of simplicity, we assume here that b and all its derivatives are bounded. 
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The motion of the fluid is described by the following system of equations 



f d t V + V ■Vx, z V=-ge z -V x ,zP 
Vx,, • V = 
Vx, xV = 



dtV - Vl+Wxvf n+ ■ = 



in fi t , i > , 
in fi t , * > , 
in fit, t > , 
for i > 0, le M d 
for t>0, X eR d 
for t > 0, X G M d 



where n 



■+ 



= (— V?7, 1) T denotes the outward normal vector to the surface and n_ 



(0.2) 



V 1 +i V6 i 



r(V6,-l) 



denotes the outward normal vector to the bottom. The first equation corresponds to the Euler equation for 
a perfect fluid under the influence of gravity (which is characterized by the term — ge z where e z denotes the 
base vector along the vertical component). The second and third one characterize the incompressibility and 
irrotationnality of the fluid. The fourth and last ones deal with the boundary conditions at the surface and the 
bottom. These are given by the usual assumption that they are both bounding surfaces, i.e. surfaces across 
which no fluid particles are transported. As far as the pressure P is concerned, we assume that it is constant at 
the surface by neglicting the surface tension. Up to a renormalization, we can assume that it is equal to zero 
at the surface. 



In this paper, we use the Bernoulli formulation of the water-waves equations. The conditions of incompress- 
ibility and irrotationnality ensure the existence of a potential flow <f> such that V — Vx.z^- From now on, we 
separate the transverse variable lei 11 and the vertical variable z € R : the operators V and A act only on 
the transverse variable X e K d so that we have V = V(f> + d z 4>. The use of the potential flow cj> instead of the 
velocity V leads to the following formulation of (|0.2[) : 
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' dt4> + \ [ | V0| 2 + \d z cp\ 2 ]+gz = -P in fi t , t > , 
A</> + <9 2 = in fi t , t > , 

d t r, - ^l+WW d n+ 4>l= v (t,x) =0 for t > 0, X £ R d 



0- 



l ii- ( P\z=-h +b(x) 



= 



for t > 0, X G 



(0.3) 



where we used the notations 9n_ = n_ 



and 9n, = n + 



V 

a, 



Separating the variables X and z in the boundary conditions and taking the trace of (|0.3|) on the free surface 
thus leads to the system : 



A(j) + d 2 ct) = 
9 t + i [|V0| 2 + |d z </>| 2 ] +g V = 
d t ri + Vri-V(j)-d z (j) = 
V6 • - 9 a <£ = 



in fi t , t > , 

at 2 = X G K d , £ > , 

at z = rj(t,X), X e R d , t > , 
at z = -/i + 6(X), X G R d , t > 



(0.4) 



We now perform a non-dimensionalisation of these equations using the following parameters : A is the typical 
wavelength, a the typical amplitude of the waves, ho the mean depth of the fluid, bo the typical amplitude 



of the bottom, to = 



/gho 



a typical period of time (y/gho corresponding to sound velocity in the fluid) and 



/gho- Introducing the following parameters : 



e= — • 0= — ■ S = — 

ho ho h 



and taking the Stokes number S to be equal to one, one gets for the non-dimensionnalized version of (|0.4 



eAcj) + d 2 (j) = 
dt<t>+\ [e|V0| 2 + |d 2 0| 2 ] +gv = 
9 t r; + eV77-V0-ia 2 = O 
d z <j)-ef3Vb- V0 = O 



-1 + j3b < z < er], X G R d , t > 

at z = £77, X G R d , t > , 

at z = £77, X G R d , t > , 

at z = -1 + /3b, X G M d , i > . 



(0.5) 



The final step consists in recovering the Zakharov formulation by reducing the previous system (|0.5|) to a 
system expressed at the free surface. To this end, we introduce the trace of the velocity potential 4> at the free 
surface, namely ip : 



^(t,X) = cf>(t,X,sr ) (t,X)) 
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and the operator Z e (eT],(3b) which maps ip to d z cf>\ z=eri . This operator is defined for any / € (C 1 n W 1 ' co )(M d ) 
by: 

/ J ffi(M d ) — > H^(R d ) \ 



Ze(e V ,Pb)f 



f 



d 



with u solution of 



V 



eAu + d^u = 0, -1 + (3b < z < er\ . 
d z u-s(3Wb-Vu = Q, z = -l + (3b 
u(X,erf) = f, X eR d . 



(0.6) 



/ 



Using this operator and computing the derivatives of ip in terms of iji and 77, the final formulation (So) of the 
water waves problem reads : 



(So) 



d t M> - ed t r)Z e (eri,l3b)i> + \ [e|V^ - e Vr}Z e (sr), (3b)i(>\ 2 + \Z £ (er,, /%)# 
dt-q + eVi]- [Vip - eVrjZ E (er], /3b)ip] = \ Z e (er), (3b)ij) . 



+ T] = 



(0.7) 



Organization of the paper 

The aim of this paper is to derive and study two different asymptotic regimes based each on a specific 
assumption on the parameter p which characterizes the topography of the bottom. The first assumption deals 
with the case (3 = 0(e) which corresponds to the physical case of a bottom with small variations in amplitude. 
The second one deals with the more complex case j3 = 0(1) which corresponds to the physical case of a bottom 
with high variations in amplitude. 

The following part will be devoted to the asymptotic expansion of the operator Z e (sr),f3b) in the two regimes 
mcntionned above. To this end, a general method is introduced and rigorously proved which aims at deriving 
asymptotic expansions of Dirichlct-Ncumann operators for a large class of elliptic problems. This result is then 
applied in each regime, wherein a formal expansion is performed and an asymptotic Boussinesq-like model of 
(|0.7|) is derived. The second and third part are both devoted to the derivation of new classes of equivalent 
systems, following the strategy developped in [6]. In the end, completely symmetric systems are obtained for 
each bottom topography regime : convergence results are proved showing that solutions of these symmetric 
asymptotic systems tend to associated solutions of the water waves problem. 



1. Asymptotic expansion of the operator Z £ (sri,f3b) 



This section is devoted to the asymptotic expansion of the operator Z e (er],(3b) defined in the previous section 
as e tends to zero, in both regimes (3 = 0(e) and (3 — 0(1). To this end, we first enounce some general 
results on elliptic equations on a strip : the final proposition gives a general rigourously justified method for 
determining an approximation of Dirichlet-Neumann operators. This result is then applied to the case of the 
operator Z e (er),/3b) and two asymptotic models with bottom effects are derived. 

1.1. Elliptic equations on a strip 

In this part, we aim at studying a general elliptic equation on a domain Q given by : 

n = {(x, z) g m. d+1 /x e m. d , -h + b(x) < z < r](x)} , 



where the functions B and r\ satisfy the following condition 
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3 h m in >0,VIeR d , r)(X) - B(X) +h > h min . (1.1) 
Let us consider the following general elliptic boundary value problem set on the domain £1 : 

-Vx ;Z .PVx z u = inn , (1.2) 

M U=, ( x)=/ and d nU\ z= _ hQ+B{x) = , (1.3) 

where P is a diagonal (d+ 1) x (d+ 1) matrix whose coefficients (pi)i<i <d+i ar e constant and strictly positive. 
Straightforwardly P is coercive. We denote by d n u \ z =-h +B(x) the outward conormal derivative associated to 
P of u at the lower boundary {z = —h + B(X)}, namely : 

d ^ u L=- h0+B(x) = - n - ' P ^x,zU lz= _ ho+mx) , 

where n_ denotes the outward normal vector to the lower boundary of f2. For the sake of simplicity, the notation 
d n will always denote the outward conormal derivative associated to the elliptic problem under consideration. 

Remark 1.1. When no confusion can be made, we denote V x by V. 

As in [6, 16, 21] we transform the boundary value problem (|1.2pp.3[) into a new boundary problem defined over 
the flat band 

s = {(x, z) e R d+1 /x e R d , -l < z < 0} . 

Let S be the following diffcomorphism mapping <S to fi : 

S : ( (X,z) .— ► s(X,z) = (r ] (X)-B(X) + ho)z + V (X)) ' (L4) 

Remark 1.2. As shown in [16], a more complex "regularizing" diffeomorphism must be used instead of (5) to 
obtain a shard dependence on rj in terms of regularity, but since the trivial diffeomorphism (<S) suffices for our 
present purpose, we use it for the sake of simplicity. 

Clearly, if v is defined over fl then v — v o S is defined over S. As a consequence, we can set an equivalent 
problem to (|1.2|)(|1.3|l on the flat band S using the following proposition (see [15] for a proof) : 

Proposition 1.3. u is solution of i i.£j )) TOj) if and only if u — uo S is solution of the boundary value problem 

-Vx, z .£VxzW = in S , (1.5) 
U =0 = / and d n u lz=1 =0 , (1.6) 



u 



where P(X, z) is given by 



P(X,z) = — i M T PM 

rj + ho — B 



W UhM(X,z)={ (V + h -B)I dxd -(z + l)V V + zVB 
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Consequently, let us consider boundary value problems belonging to the class (|1.5[)(|1.6[) . From now on, all 
references to the problem set on S will be labelled with an underscore. 

On the class (ll.5p(|1.6p of problems set on the flat band S, we have the following classical existence theorem : 
assuming that P_ and all its derivatives are bounded on S, if / G H k+ z (R d ) then there exists a unique solution 
u £ H k+2 (<S) to (11. 5|) (|1.6j) . The proof is very classical and we omit it here. 

As previously seen, we need to consider the following operator Z(rj, B) which maps the value of u at the upper 
bound to the value of d z u\ z=v : 



F§(R d ) — ► H?(R d ) 

f i — ► d z u\ z=v with u solution of (| 1 . 2[) (11.3 



Z{r),B) 



Remark 1.4. The operator Z £ defined in (J076J) corresponds to the operator Z in the case where P 
in (0)01). 



eh 
1 



To construct an approximation of this operator Z{rj^ B), we need the following lemma which gives a coercitivity 
result taking into account the anisotropy of (|1.2|)(|1.3|) . 



Lemma 1.5. Let r, € W 1 ' 00 (R d ) and B € Fy 1 '°°(R <i ). Then for all V € R d+1 : 

(V,PV) > co(||»?|ki.-,||fl|ki.=o) \Vpv\ 2 , 
where Co is a strictly positive function given by 

. Pd+i 

i I i mm 

/ \ 1mm i 1 l<i<d pj 



(d + l) 2 h min (x + h + y) {x + y)< 



Proof. Using Proposition 1 1.31 . we can write, with S(X) = r/(X) + ho — B(X) 

1 



(V,PV) = 



V , M 1 P MV 
MV,PMV 



) 



i VP MV,VP MV 

±m(Vpv)\ 2 

where M. — \fP M(^/P~)~ 1 . Thanks to the condition we deduce the invertibility of M and hence the 

invertibility of M.. This yields the following norm inequality for all U S R d+1 : 



|LT|<(d+l) 



VdM- 



-=MU 
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with 



where \A\ C 



M.- 1 = 



1 T 





i 



sup \ a i.j\L°^(R d ) an< i P ls the d x d diagonal matrix whose coefficients are {pi)\<i<d- 

\<i,j<d+l 

If we apply the previous inequality to our problem, one gets : 



(V,PV)> 



1 



VSM' 



P V 



(d + 1) 2 

Thanks to the expression of .M -1 given above, we obtain the following inequality : 

(V , PV) > co(\\v\\w^,\\B\\ w ^) \VPV\ 2 , 
where cq as in the statement of the Lemma 11.51 



□ 



Let us introduce the space H k '°(S) : 

H k <°(S) = {vEL 2 (S), \\v\\ h h.o := (J W-,z)\g hm dX <+co} 



The result of this subsection consists in the following theorem which aims at giving a rigourous method for 
deriving an asymptotic development of Z(r], B). Of course, P, and thus P, as well as the boundaries r\ and B, 
can depend on e in the following theorem. In such cases, the proof can be easily adapted just by remembering 
that < e < 1. 



Theorem 1.6. Let p e N*, k e W, T] e W k+2 >°°{R d ) and B G W k+2 '°° (R d ) . Let < e < 1 and u app be such 
that 

-V x ,z-PV x , z u app = e p R e in S , (1.7) 

t^app \ z= q / > 9 n U app \ z=1 — £ P T , (1-8) 



where (R £ ) o<e<i an d ( r£ )a<e<i ar ^ bounded independently of e respectively in H k+1 '°(S) and H k+1 (M. d ) 
Assuming that h m i n is independent of e and that tht 
bounded by a constant 7 independent of e, we have 



Assuming that h m i n is independent of e and that the coefficients (pi)i<i<d+i of P are such that )i<i<d are 



Z(v,B)f 



1 



rj + h Q - B 



(d z u app ) | a!=0 



< 



H k +h y/Pd+i 



C fe+2 (||i? E |U +1 ,o + |r e |^ +1 ) 



where Ck+2 = C(\r]\yyk+2,ao , |i?|jyk+3,«>) and C is a non decreasing function of its arguments, independent of the 
coefficients [pi)i<i<d+\- 
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Proof. In this proof, we often use the notation Ck — C(\i]\ W k,oa , |S|^fc,»o , ho, h m i n , k, d, 7) where C is an unde- 
fined non decreasing function of its arguments. The notation Ck can thus refer to different constants, but of 
the same kind. 

A simple calculus shows that Z(rj, B) can be expressed in terms of the solution u of (|1.5[)(|1.6[) via the following 
relation : 

Z(r>, B)f = -i d z m 

W ' u r, + h Q -B ~\*=° 

Using this fact, we can write 

Z(7? ' B)f r 1 + h -B d * U °^o = v + l_ B ^ I- ■ 
Introducing <p :— u app — u we use a trace theorem (see Metivier [20] p. 23-27) to get 
Tif- 

77 + h Q - B 



Z( V ,B)f - 1 (d z u app ) U J Hk+ i < C k+1 (\\d z <p\\ Hk+1 , + \\d 2 z <p\\ Hk , ) . (1.9) 



It is clear that the proof relies on finding an adequate control of ||9 z y>||#*+i,o and ||d|y>||#fe,o. The rest of this 
proof will hence be devoted to the estimate of both terms. 

1. Let us begin with the estimate of ||9 z (^||jyfe+i,o. To deal correctly with this problem, we introduce the 
following norm ||*||rji defined by : 

IMIffi : = IIVpVx.^IU 2 ^) • 

First remark that for all a € N d such that |a| < k, d a ip solves : 

f -Vx,z-PVx.zd a <p = e*'d a R e + V x ,z-[d a ,P]V x ,zV , 

{ pi (1-10) 

{ ff*<p\ t=0 = 0, ]M= - 1 +d [B n ' Ei <p\ t= -i=e*d a T* . 

In order to get an adequate control of the norm ||i9 z (/3|| ff fc+i,o, we prove the following estimate by induction on 
\a\ < k: 

VaeN d /\a\<k, \\d a ^\\ H1 <-^C k+1 {\\R e \\ Hk ,o + \r"\ Hk ) . (1.11) 
The proof of (|l.lip is hence divided into two parts : initialization of the induction and heredity. 
• Initialization : |a| = . 

Taking a = 0, multiplying (|1.10|) by ip and integrating by parts leads to : 



(S) 



The boundary term at the free surface vanish because of the condition <p\ z= o = and using the condition at 
the bottom leads to : 

(£Vjc,*¥>, V x ,z<p) L 2, s) = (e p R £ , v) L 2 {s) +e p [ rV|. = _! ■ 

Finally, using Cauchy-Schwartz inequality, one gets : 

(£Vx, z V, Vx^) i2(s) < e p ||i? e ||L2(5) +£ p Ir^l^^d) |^ u= _ 1 | i 2 (K d ) . (1.12) 
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Recalling that **p\ z=a — and that the band S is bounded in the vertical direction, one can use Poincare 
inequality so that < H^Hl^s) and \ip\ z= _ 1 \ L 2 {Rd) < \\d z ip\\ L 2 {S) . Therefore, (|1.12|) yields 

(PV x ,z<P,Vx,zV) LHs) < -^=\\R e \\ms)\Mm + ^=\r £ \mm«)\MHi ■ (1-13) 

yPd+1 yPd+i 
Using Lemma ll.5l to bound (£Vx jZ f, Vx i2 <p ) L 2/g\ from below, one finally gets : 

co(Hw^,\B\ w ^)M\l < (\\R e \\H°-° + \r e \ H o) |M| ffl . 

\/Pd+l 

Since co( |ry|v7i,«>, l-Blvi/ 1 . 00 ) depends only on h m i ni d and 7 through the quantity nrmi<j<d (by Lemma ll.5|) . 
and since the function cq is a decreasing function of its arguments (again by Lemma 1 1 . 5 j) , we get the following 
desired estimate : 



p 



IMI^ < -£=Ci (\\R e \\ H o,o + \r e \ H o) 

\/Pd+l 

which ends the initialization of the induction. 

• Heredity : for m £ N* fixed such that m < k, we suppose that (|1.11[) is verified for all a £ N d such that 
\a\ < m — 1. 

Let a £ N d such that \a\ = m. Multiplying (11.10|) by d a ip and integrating by parts on S leads to : 



(5) 



The boundary terms at z = vanish because of the condition d a (p\ z= o — 0, and using the second boundary 
condition d n (d a (p) i z= _i + 9 „ '^"V|*=-l = £ p d a r e , one gets : 

and with Cauchy-Schwartz : 

(£Vx,,0°V, Vx,,9» i3(s) < ^||fl a ie 6 ||^(s)||5 vlU»(s)+e p |fl a r B | La(Rd) |a a ^|, = _ 1 | £ »(R-) 

+ |([5 Q ,£]Vx,^, Vx, z 5>) i2(s) . 

By using the same method and arguments as in the initialization, the following inequality arises : 

co( H w i.», \B\ wl , x )\\d a <p\\l < -j?= (\\RZ\\ Hk ,o + \rt\ Hk ) ||9>||^ + |([9 a ,P]Vx,^, V x , 2 fl a¥ > ) l2(5) 

(1.14) 

Let us now focus on the second term of the left hand side of p. 141) . In order to get an adequate control of this 
term, we have to write explicitly the commutator [9 Q ,P] : 



[d a ,P]V x ,z<P = C(\a'l\a"\)d a 'PV x , z d a " 



a. +a —a 
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where C is a constant depending only on \a'\ and \a"\. This leads to the expression 

(F,Z]Vx^, Vx,9» i2{5r C(\a'\,\a"\)(d a 'PV x , z d a "ip, V x , z d° 



a -\-Oi —a 



L 2 (S) 



From now on, we just consider a single term of this sum. Using Proposition 1 1 . 31 we derive the explicit expression 
of P and deduce from it the explicit expression of d a P_ : 



d a P = 



d a t? - d a Bj P d P d d a U 

P d d a 'u) T d 



T 

a ( Pd+l+U-Pd U 



ri+ha-1 



where P d is the diagonal (d x d) matrix whose coefhcents are (pi)i<i<d, and U the vector defined by U = 
— (z + l)Vr/ + zVB. Using this expression, one easily gets (with V = Vx) : 



d a PV x ,zd a <p, Vx ^ da <P) L2[s} = [{d a V- d a B)P d Vd a tp , Vd a <p) ^ ^ 
d z d a "ipP d d a 'U , Vd a w) +( P d d a 'U ■ V<9 Q 'V , d z d a ip) 

> L 2 (S) \ ) L 2 (S) 

If we focus on the first term of the right hand side of this equality, we easily get the following intermediate 
control using Cauchy-Schwaxtz inequality and the definition of ||*||rji • 



((d a ' v -d a 'B)p d vd a "v, va>) < + \B\ wla >^)\\V^vd Q ''v\\ LHS) \\y/F d vd Q v\\ LHS) 

\ / L 2 (S) 



'L 2 (S) 

< (\ri\ w *,~ + \B\ w >.ao)\\d°''<p\\ jil \\& a <p\\ lil 



< 



C k+1 (\\R £ \\m.o + \r^\ Hk )\\d a v \\ Hl 



To derive the last inequality, we used the induction hypothesis on ||(9 Q <p\\r}i since \a"\ < m — 1. 

Let us now focus on the second term of the right hand side of (|1.15p . Using the same arguments as previously 

and Poincare inequality, the following controls arise : 

(d z d a " v p d d a 'u, vd>) < \\^ d d a 'u\\ 00 \\d z d a ''^\\ L 2 (S) \\^p~ d vd a ^\\ L 2 [s) 



< 



Pd+l 



< -^=c k+1 (\\ir\\ H H. + \t*\ H u)\\d a <p\\ lil , 



since - Pd - °° < 7. 

Pd+i — 1 

The control of the third term of the right hand side of (|1 . 1 5f> comes in the same way : 

(p d d a 'U-Vd a "v, d z d a y) <^ = C k+1 (11^11^.0 + 1^*) . 

\ / L 2 (S) ^/Pd+1 
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The next step consists in controling the last term of the right hand side of (|1.15[) . We need to do some preliminary 
work on this term before attempting to estimate it adequately. A straightforward calculus gives us : 

A/o 

+ U-P d Ud a ' (—± -)+p d +id a,/ ' 



v 7] + ho — B ) \ r\ + ho — B 

We plug the previous writing and use the same tools as previously to get : 

(V ( Pd+1 *?' P Z U ) d z d a "v, d z d a v) < \\P d \\ooC k+1 \\d z d a " y\\ L * {s) \\d z d a ^\\ L 2 (s) 

V V V + ha-B J J L2(s) 



+ C k+1 \ \y/p d +id z d a <p\\ L '(s) \ WPd+id z d a ip\\ L 2 (s) 

Mkc, +1 ||a Q 'V||^||9>||^ 

Pd+l 

+C k+1 \\d a "v\\ lil \\d a v\\m 



< 



p 



^Pd+i 



C k+1 (\\R £ \\m-o + \rZ\ Hk ) na^n^ ( 



where we once more used the induction hypothesis. 

Gathering the four previous estimations of each term of the right hand side of (I1.15|) and using the explicit 
writing of the commutator [d a ,P] leads to the final estimate of |([<9 a ,P] \?x, z f , Vx.z c*" 1 / 5 )l 2 (S)| : 

([d a ,P]Vx, z <fi, V x ,*d<* V ) L2{s) \<^=C k+1 (\\mU.o + \^\ Hk ) . 
The last step simply consists in pluging this last estimation in the estimation (|1.14[) . which gives : 

co(\v\w^,\B\wi.-)\\d a <p\\l <-^||i? £ |Uo||9>||^ . (1.16) 

■s/Pd+l 

As in the initialization, this last estimation leads to the desired result, which ends the heredity and hence the 
proof of (TTTTj) . 

To conclude this first part of the proof, we use the fact that : 

\\d z <p\\ H k+i,o < C(k + 1) sup \\d z d a <p\\ L 2 {s) 

\a\<k+l 

C[k +1) ina || 

< . SUp \\d ipWjji , 

y/Pd+1 \a\<k+l 

and the estimate we just proved to finally get : 

HS^II^+i.o < — C k+2 (\\R E \\ H * + i,o + \r"\ Hk+1 ) , (1.17) 
Pd+i 

which ends the first part of the proof. 



TITLE WILL BE SET BY THE PUBLISHER 



13 



2. In this second part, we aim at controlling the quantity | \d^tp\ |ifk,o . To this end, we prove with a direct 
method the following estimate : 



\\d 2 MU.o < C k+2 (\m\ H *+r,o+\r s \ Hk+ r) 

Pd+l 



(1.18) 



We first use the equation satisfied by ip in order to express d z (p in terms of other derivatives of ip such as 
V(/?, d z Lp, Vd z (p or Aip. There comes the following expression : 



,2 / V + ho-B \ . _ „ d z {U-P d U) 



Kv = 



Pd+1 +U-P d U 



[-Vx,, ■ Q Vx^ - z ' h a ' d z <p ~ # R s \ , 
r/ + ho — B 



, n ( (v + ho-B)P d P d U 
The following estimates arise (using ||u||jy*,o < C{k) sup| a | <fe ||9 a u||i2( lS ) ) : 



\\d 2 M\ m , < C k \\^x,z-Q^xM\HKo+C M \\P d \\ 00 \\d t ip\\ H H.o+e p \\B! 1 \\ H ufi 

Pd+i 



< Ck+i 

Pd+i 



C{k) sup \\d a {V x , z -QV x , zV )\\ L 2 (s) 

|a|<fe y/Pd+1 \a\<k 



C(k) sup 



\H k .° 



< — C k +i sup \\d a (Vx.z ■ QV X ,^) \ \l'(S) + —— -L^^ C k +i {\\R e \\ H k, a + \r £ \ H k) 

Pd+l \ a \<k Pd+1 Pd+1 



Pd+1 
1 



Ck+i \\R 



+1 1 1 -ft fffc.o 



< Ck+i sup \\d a (V x , z -QV x ^) \\ L '{S) + C k +i(\\R e \\H^ + \r s U) , 

Pd+l \a\<k Pd+l 



where we used the result (fTTTTl) and the fact that l|fj||o ° < 7 



Pd+l 



The last part of the initialization aims at correctly estimating the norm \\d a (Vx,* ' Q^x,z^p) ||l 2 (s)- The 
explicit writing 



(V x>z -QV x , z <p) = C{\a'l\a"\) (y x .z-d a 'Q^ x , z d a " 



'•P 



a.' -\-a f/ — a. 



and the expression of Q furnishes us with the following estimates 



14 



TITLE WILL BE SET BY THE PUBLISHER 



\\v x , z -d a Qv x , z d a cp\\ LHS) < c W \ \\v-p d vd a <p\\ L 2 (s) 

+c [a > l+ i ||V^II«>||V^w*'V|| i2(s) 
+c M+a II^IUH^'Viu^s) 
+q Q 'i+i ||P<illoo||^V9 a 'V|U2 (s) ) 

< c k+2 ( ||v-p d va Q 'Vll^(s) 

+||V^||oo l|9 a 'V||^ + ^feirVll^ ) 

< C k+2 (||V-P d V9 a 'VlU=(s)+ ^Cfe+i (||ii E ||^,o + |r e |^)) 

We estimate the term ||V • Pd\7d a ip\\i, 2 (S) using the following technique : 

||v-p d va«Vll^ ( 5) < llV^IU £ llvKC^'Vll™ 

l<t<d 

< iiv^iu £ ll^'VlU 



< dllV^Hoo sup ||0"V|| ril 

|m| = |a"| + l 



(1.19) 



< d||VPd||oo sup ||3 m ¥ >||^ 1 

|ro|<fc+l 

< ePC k+2 (II^II^+lo + \r"\ Hk+ i) . 

We plug this result in (|1.19p to obtain 

\\Vx,z-d a 'QVx, z d a "^\\ L 2 (s) <e p C k+2 (\\R E \\ Hk+ i,a + \r s \ Hk+1 ) , 

which finally leads to 

\\d a (Vx,z-QV x , zV )\\ms) <e p C k+2 (II^II^+lo + I^Ih^+O ■ 
This way, we get our desired estimation of ||S|</?||#fe,o: 

H^IIhw < —C k+2 (1^11^+1.0 + |r £ | ff , +1 ) . 
Pd+i 

Gathering (fTTT7|) and (fTTTH)) in (fl~9]) ends the proof of the theorem. □ 
1.2. Application 

We recall that by definition, Z s (er], (3b) f = d z u\ t= where u is solution of the boundary value problem 

eAu + d*u = inn , (1.20) 
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«!.«,=/, (»»«-e/3V6-V«) u= _ 1+/Jb = , (1.21) 

This elliptic problem p.20|) (|1.21[) belongs to the class of general elliptic problems (|1.2p (|1.3[) defined in the pre- 
vious subsection. The corresponding matrix P is here designed by P e : 

P £ = ( £ld Q xd J ^ . (1.22) 

The upper boundary of f2 is here defined by {z = er/} and the lower one by {z = — 1 + fib}. We make the 
additionnal assumption that e and are bounded in the following sense : < e < 1 and there exists a strictly 
positive constant (3q such that and < j3 < (3q. Furthermore, condition Ijl.ip is here verified thanks to condition 
(jO.ip . And finally, we remark that ) 1<i<d are bounded by 1 since < e < 1. Our goal is here to apply the 

Pd-\-l 

previous theorem to get asymptotic estimates on Z e (er],f3b). 

We recall that we are here interested in two differerent regimes depending on the (3 parameter. The first one, 
namely (3 — O(e), refers to the physical case of a bottom with variations of slow amplitude. The second one, 
namely (3 = O(l), refers on the contrary to variations of high amplitude of the bottom. In order to improve the 
readability, we take /3q = 1 : we thus write (3 = e for the first regime and (3=1 for the second one. 

1.2.1. The regime (3 — e : small variations of bottom topography 

The boundaries of the domain Q are here defined by {z = er]} and {z = — 1 + eb} while matrix P £ remains 
as in (|1.22[) . Thanks to Proposition 11.31 we are able to set an equivalent problem to p.20[ip.21[) defined over 
the flat band S : u = uoS then solves the problem : 

- V X:Z • P e Vx,z u = in S , (1.23) 
% =0 =/ , 9 nMu= _ 1= . (1.24) 

where the matrix P_ £ is given by 

e(l + £(r,-b))I dxd -e 2 [(z + l)Vry - zWb] 

The following result gives a rigourously justified asymptotic expansion of Z £ (erj,(3b)f as e goes to : 

Proposition 1.7. Let k e N, r? € T^ fc+ ^°°(]R d ) and 6 e VF fc+2 '°°(M d ). 
T/ien /or / such that V/ G H k+6 (R d ), we have : 

\Z e (en,(3b)f-(eZ 1 +e 2 Z 2 )\ Hk+1/2 < e 3 C k+2 |V/|^ +6 , 

with : 

f := -A/ , 

\ Z 2 := -|A 2 / - (77 - 6)A/ + V6 ■ V/ . 



P e = 
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Proof. To prove this proposition, we use essentially Theorem 11.61 with p = 3. We know that ( are 

\ Pd + 1 /i<i<d 

bounded by 1. Thus, in order to derive an asymptotic expansion of Z £ (eri,(3b)f, we only need to compute an 
approximate solution u app which satisfies the hypothesis of Theorem 11.61 for p = 3. This approximate solution 



form M a; 
to say : 



can be constructed as in [6] using a classical WKB method, which consists in looking for u app under the 
eu± + e 2 U2- We want this function to verify the properties required by Theorem II. 6[ that is 



u Q 



-Vx,z-P £ Vx, z u app = e p R £ in S 

W'app \z = a f i ^app\ z = ^i £^T 



(1.25) 
(1.26) 



where (i? e )o< e <i and (r £ )o <e< i are bounded independently of e respectively in H +1,0 (S) and H k+1 (M. d ). 
We decompose the matrix Pf_ under the form = Pq + eP\ + e 2 P2 + e 3 P e where Pq, Pi, Pi are independent 
of e, and if we plug the desired expression of u app into this problem, we get R £ = V • T £ and r e = e z • T £ _ 
where T £ = P^S x,zU\ + Pi ^x,zU2 + P £ Vx z (m + u\ +1*2), and the following system of equations and boundary 
conditions on uq, u\, u-i : 

d 2 u = , 

9>i + (A - (77 - b)dfj u = , 

9> 2 + (A - ( v - b)df) u x + (rj - b)Au -2[(z + 1)V/ - zVb] ■ 

vd z u - [{z + i)A/ - zAb] ■ d z u -( n - b) 2 d 2 z u = , 



u ou =0 = / 



with < 



l i\z=o 



= 0, 1< i < 2 



d z u iU= _ x =0, < i < 1 , 
. 9 z u 2 \ z= _ 1 - V6 • Vw | i= _ 1 = . 

We can verify that the following values of 1*0,1*1, it 2 satisfy the previous equations and boundary conditions : 



wo = / , 
1 
2 



141 



"2 



Z + l) 



A/ 



Using these values of uo,i*i and 1*2, one can easily obtain the following estimates of R £ and r £ : 

\\R e \\h^> < C k+ 2 \Vf\ Hk+e , \r £ \ Hk+1 < C k+2 |V/| Hfc+3 . 



Thus 



1 o pp 



satisfies the properties required to apply Theorem 11.61 The last steps of the proof consists in com- 



puting 1+£ ^ v _ b j dz u app\ z=Q using the explicit expression of u app previously determined, and then apply Theorem 
11.61 An easy Taylor expansion of 1+e ^ v _ b ^ d z u app \ z=0 then yields the result. □ 



Remark 1.8. The method developped here to get and prove our asymptotic expansion is improved compared 
to the one developped in BCL since we do not need here to compute the term 1*3. 



TITLE WILL BE SET BY THE PUBLISHER 



17 



Remark 1.9. If we take b — - i.e. if we consider a flat bottom - , we of course get the same expansion as the 
ones proved in [6]. 



I. 2.2. The regime (3=1: strong variations of bottom topography 

The boundaries of the domain O are here defined by {z = erj} and {z = —1 + b}. Using again Proposition 

II. 31 we set an equivalent problem to (|1.20|)(|1.21|) defined over the flat band S : this new problem is the same 
as the one defined in the first regime, at the exception of the matrix Pf which is now given by 



P £ = 



e(l + en-b)I dxd -e[e(z + 1)Vt? - zVb] 
-e[e(z + l)Vi7 - zVbf i+sMs+^YMl 

As in the first regime we give a rigourously justified asymptotic expansion of Z e (er], (3b) f in the present regime. 

Proposition 1.10. Let k e N, 77 e TU ,£+2 ' 00 (R' i ) and b G W rfc+2 ' 00 (IR d ). 
Then for all f such that V/ S H k+6 (R d ), we have : 

\Z £ (er 1 ,(3b)f-(eZ 1 +e 2 Z 2 )\ Hk+1/2 < e 3 C k+2 |V/|«»+« , 

with : 

' Zx :=-V- ((l-b)V/) , 

Z 2 := i V • - bf VA/ - (1 - bf VV • ((1 - b) V/) ) - nAf . 



Proof. The proof of this proposition follows exactly the same steps as the proof of Proposition ll.7l The following 
values of uq, Ui, u 2 are found : 



U = 


f 




Ul = 


(1 






2 




(1 


-bf 


u 2 = 


24 



(l-(z + lY)Af + z(l-b)Vb-Vf , 
A 2 / z 4 + H-l-L AV • ((1 - b)Vf) z 3 - (1 - b)nAf z 2 



6 

"^y^ V • (^^VA/ - (1 - 6) 2 VV • ((1 - b)Wf 
-77(2(1 -6)A/ + V6- V/) z . 

The error bound is derived in the same way and the previous values lead to the result. □ 



Remark 1.11. By formally taking b = eb, we recover the result of Proposition 1 1.71 
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1.3. Derivation of Boussinesq-like models for uneven bottoms 

We recall the Zakharov formulation of the water waves equations, from which we intend to derive new systems 
using results of the previous subsection : 



(So) 



d t il>-ed t r)Z e (er),/3b)i> + ± e \ Vip - e V 7 /Z e (er/, /3b)^\ 2 + \Z e (en, (3b)^\ 2 + V = , 
dtri + eVrj- W - eVr)Z e (eT), pb)ip = - Z £ (en, (3b)ip . 



As in [6] , we introduce the notion of consistency. 

Definition 1.12. Let a, s £ K, e > 0, T > and let (V s ,if) 0<e<eo be bounded in W l <°°([0, f ]; H a (R d ) d+1 ) 
independently of e. This family is called consistent (at order s) with a system (S) if it is solution of (S) with a 
residual of order e 2 in L°°([0, §]; H s (R d ) d+1 ). 

Thanks to the previous results, we are now able to enounce the following propositions which show the consistency 
of two Boussinesq-like systems with the system (So). We introduce here the quantity h = 1 — 6 which corresponds 
to the non-dimensional still water depth. From now on, this quantity is considered as a bottom term since it 
only depends on b. 



Proposition 1.13 (Small variations regime j3 = e). Let T > 0, s > 0, a > s and (ip £ ,i] £ )o<e<Ea be a family 
of solutions of |0.7| j such that (Vip 6 , f? e )o<e<£ * s bounded with respect to e in H /1 ' oo ([0, — ]; H a (R d ) d+1 ) with a 
large enough. We define V s := \7?p 6 . Then the family (V s ,rj)o<e<e * s consistent with the following system : 



(Bi) 



d t V + V77+ §V|T/| 2 = , 

dtv + V • V + e [v • ((rj - b) v) + ± AV • V 



= 



Proof. This is clear thanks to the asymptotic expansion of the operator Z e (er},(3b) : plugging this in system 
(|0.7|) . neglecting the terms of order 0(e 2 ), and taking the gradient yields the result. □ 



Proposition 1.14 (Strong variations regime (3 = 1). Let T > 0, s > 0, a > s and (ip £ , i] e )o<e<eo be a family 
of solutions of |fl.7| j such that (VV> £ , f? £ )o<e<£ * s bounded with respect to e in W 1,o °(\0,—]; H a (R d ) d+1 ) with 
a large enough. We define V £ :— Vip E . Then the family (V s , ?y)o<e<e * s consistent with the following system 
( with h = 1 — b) : 



d t V + V?7 + §V|^| 2 = 



(Bi) 



dtV + V • (hV) + e V • (ijV) - iV • VV • V - h 2 VV ■ (h V) 



= 



These results close the first part. From now on, our work is divided into two parts, each corresponding to one 
of our two regimes f3 — e and (3 = 1. Each analysis starts from the two previous Boussinesq-like models. 
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2. The regime of small topography variations 



We recall our previously derived Boussinesq-like system (B\) on which we base our analysis : 

' d t v + vi 1 + ^v\v\ 2 = o , 



(Bi){ 



dtv + V • V + e [V • ((r? - b) v) + \ AV • V 



= 



We now follow the method put forward in [5] and [6] to derive equivalent systems to (B\) in the meaning of 
consistency. The rigourous justifications of the derivation of these systems is adressed in section 2 in [6]. 

2.1. A first class of equivalent systems 

As in [5] (for the ID case) and in [6], we define : 

Ve= (l + |(l-e 2 )A) V , 

which corresponds to the approximation at the order e 2 of the horizontal component of the velocity field at 
height -1 + 9 for 6 E [0,1]. If we remark that Vg = (l + §(6» 2 - 1)A) _1 V + 0{e 2 ), the expression of Vg in 
terms of V comes in the following way by supposing V regular enough : 

V= (i+ £ -(9 2 -1)a) V e + 0(e 2 ) , 

where 0(e 2 ) is to be taken in the L°°([0, ^},H s (R d )) norm. 
Plugging this relation into the system {B\) leads to : 

d t V + V V + ^(v\V\ 2 + (6 2 - l)Ad t Vo^j = 0{e 2 ) , 
dtv + V • V e + e [v • ((77 - b) v) + - i) AV • V e ] = O(e) 2 . 

At this point we use the classical BBM trick which consists in writing the following approximations at order 
0(1) , coming from the previous equations : 

d t V e = -Vf? + O(s) = (1 - n)d t V B - + 0(e) , 

V • V e = -dtv + 0(e) = AV - Vg — (1 — \)d t r) + 0(e) , 

where A and \x are two real parameters. 

We plug these relations into the dispersive terms of the last system to get : 

' d t Vg + Vry + § [V|^| 2 - fi(6 2 - l)AVr? - (/* - l)(6 2 - l)Ad t V e ] - 0(e 2 ) , 

dtv + V • Vg + e [V • ((t? - b)Vg) + X (£ - i) AV • Vg - (1 - A) (£ - i) A^r?] = 0(e 2 ) . 

We then introduce the class 5 of all the systems of the previous form : these systems are denoted by Sg^ ^ and 
can be rewritten in the compact form : 



' (l-ea 2 A)d t V e + V7 1 + e [\V\Vg\ 2 + a^V-q] =0 , 
(1 - ea 4 A)dtv + V • Vg + e [v • (fa - b) Vg) + a 3 AW ■ Vg 



. 
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with 

«3 = A(f-i), a 4 = (l-A)(f-i) . 

On this class 5, the previous computations give us the following two results of consistency : 

Proposition 2.1. Let 9 £ [0,1] and (ip £ , r] £ )o <e<eo be a family of solutions of 9f ) such that (Vip £ ,rj £ )o <e<eo 
is bounded with respect to e in W l '°°([Q,^]; H a (R d ) d+1 ) with a large enough. We define V £ = V-0 £ and 
Vg = (l + |(1 — (9 2 )A) V 6 . Then for all (A,/x) € K 2 , the family (V g £ , r) £ )o <E<eo is consistent with the system 

(^0,A,/i)- 

Proof. We saw in the previous section that if (ip £ , r) £ )o <£<eo is a family of solutions of (|0.7p . then the family 
(V-0 e , V £ )o<e<eo is consistent with the system (Si). Thanks to the previous computations, and since the choice 
of the parameters (A, n) is totally free, it is clear that (V g £ , r] £ )o <s<EQ is consistent with any system (S g I ). □ 



Proposition 2.2. C/p to a change of variables, all the systems belonging to the class S are equivalent in the 
meaning of consistency. 

Proof. Let (0, A, /x) € [0,1] x R 2 and {V e £ ,n £ ) 0<e<Eo a family consistent with (S e ). We then define, for 

0i e [o,i], 

v£ - (i + |(i - <?)a) (i - |(i - 2 )a) ; 

using the fact that (l — |(1 — 2 )A) V g £ = V £ + 0(e 2 ) and the previous proposition, we easily deduce that the 
family (V^ , 7 7 e ) 0<£<£ ^ s consistent with any system (S 6i 1 Xi pi ) for any (Ai,/ii) G K 2 . □ 



Remarks 2.3. 

• By taking 9 = 1,A = 1,/i = 0, we remark that the previously derived Boussinesq-like system B± is 
actually a member of the class S 

• By taking A = fi = 1/2 and 2 = 2/3, we get a\ = a 2 = = 04 = so that the dispersive part of the 
correponding system (S g ^ ) is symmetric. However, the nonlinear terms, that are not affected by the 
choice of 9, A, fi, are not symmetric : this problem is adressed in the next section. 

• In [8], Chen formally studied in ID the case of slowly variating bottoms and derived the same class 
of systems at the exception that she considered time- dependent bottoms : her systems thus contain 
additionnal time derivative terms on the bottom that does not appear here but could be easily obtained 
for a time dependent bottom. 



2.2. A second class of equivalent systems 

Adapting the nonlinear change of variables of [6] to the present case of varying depth, we introduce V : 

V = (l + I (V - bj) V . 
This nonlinear change of variable symetrizes the nonlinear part of the equations. 
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This change of variables only affects the nonlinear terms and not the dispersive terms. If (V £ , r) e )o <e<eo is 
consistent with a system (S s ^ ) of the class S, then V £ = (1 + § if) V e and rf satisfy the following equations : 



(1 -ea 2 A)d t V £ + Vrf + e ±V|7f | 2 + \^\V e \ 2 + \V e V • V £ - \b Vrf + aiAVrf = 0(e 2 ) , 
(1 - ea 4 A) d t rf + V-V E +e [|V • ((?f - 6) f e ) + a 3 AV • V e ] = 0(e 2 ) . 



As observed in [6], if we consider a two-dimensional domain, that is to say d = 1, the nonlinear terms are 
actually symmetric. But this is not the case in a three-dimensional domain. However we can deal with this 
problem for d — 2 using the following remark coming from [15] : 

^V\V £ f = ^7\V £ f + ^(V £ ■ V)V £ + ^V £ A (V x V £ ) 

Assuming that W x V s ~ 0(e), one formally derives the following system : 



(1- ea 2 A)d t V £ + V?f + e \V\v £ \ 2 + |V|f £ | + \(V £ ■ V)V e 

+§ V £ V • V £ - \ b Vr/ e + aiAVrf 
(1 - ea 4 A) 9 t 7y e + V • V s + e ±V • ((rf - 6) V e ) + a 3 AV • U £ 



= 0( £ 2 ) , 
0(£ 2 ) • 



The nonlinear terms of the previous system are now symmetric regardless of the dimension. This previous 
computations are summed up in the following proposition : 



Proposition 2.4. Let (V s , ri £ )o <e<eo be a family consistent with a system (S e \ ) and V £ — (l + | r/ £ ) V £ . If 
V x V £ = 0(e), then the family (V £ , rj £ )o <e<£o *s consistent with the following system : 



(1 - sa 2 A) 8 t V + Vn + e ±VM 2 + ±V|U| 2 + \(V ■ V)V + i V " V ■ V - ± 6 Vt? + oiAVt? 



(1 - ea 4 A) d t r] + V • V + e ±V • ( (77 - 6) V ) + a 3 AV ■ V 



= 



We introduce the class T composed with the systems of the form (T g ^ ) for any (0, A, /i) € [0, 1] xl 2 . Using 
this result, we prove the following proposition : 



Proposition 2.5. Let 6 € [0,1] and (tp £ 7 rj £ )o <£<£o be a family of solutions of |ft 9j ) such that (Vip E , r) £ )o< £ <e 
is bounded with respect to e in VT 1:OO ([0, ^]; H a (M. d ) d+1 ) with a large enough. 

We define V £ = (l + §(ry - 6)) (l + f (1 - 6> 2 )A) V^ £ . TTien /or aZJ (A,/i) S R 2 , tte /ami??/ (V £ , ?f )o<e< £0 
consistent with the system (T g 

Proof. Thanks to Proposition ^. 1[ the family (V/, ?y e )o<e<e ^ s consistent with the system (S e | ) for all (A, /1) e 
K 2 , where VJj £ = (l + |(1 — # 2 )A) V £ and V e = Vip £ - We then use the following remark : by hypothesis, the 
velocity field V is irrotationnal, thus V £ = 0(e) and V £ — 0(e). Applying the previous proposition yields the 
result. □ 
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2.3. A new class of completely symmetric systems 

We remarked in the first section that there exists values of (8, A, /i) such that the dispersive terms are 
symmetric. Consequently, the corresponding system (Tg ^ ) of the class T is completely symmetric since both 
its dispersive terms and nonlinear terms are symmetric. We thus introduce the non-empty subclass of T denoted 
by T, composed with the systems of the form (T s £ ) for which we have a\ = 03, 02 > 0, 0,4 > 0. The first 
condition a± = 03 symetrizes the nonlinear terms and the last ones 02 > 0,04 > ensure the well-posedness 
of these completely symmetric systems. Indeed, one of the great advantages of these systems belonging to £ is 
that we have a well-posedness over a long time scale : 

Proposition 2.6. Let s > 5 + 1 and (9, A, /i) be such that the system (T e ^ ) belongs to the class E. 
Then for all (Vo,?7o) 6 H s (M. d ) d+1 , there exists a time Tq > independent of e and a unique solution (V,rj) G 
C([0,f };H s (R d ) d+1 )D C 1 ^,^- ];tf s - 3 (M d ) d+1 ) to the system (T g JJ such that (V, V ) U=0 = (V&.ifc). 
Furthermore, this unique solution is bounded independently of e in the following sense : there exists a constant 
Co independent of e such that for all k verifying s — 3fe > i + 1, we have : 



Proof. This theorem is a very classical result on hyperbolic symmetric quasilinear systems, and we omit the 
proof here. □ 



As in [6], we are now able to rigorously construct approximate solutions to the water waves problem from the 
solutions of any of these symmetric systems. 

More precisely, let us consider a solution (tp £ ,rj £ ) to the initial system (|0.7p with initial data (V'SS^o) sucn that 
(Vt/iq,^) £ H s {R d ) d+1 for a suitably large value of s. We define V s = V-0 e and Vq = Vt/'o- From this solution 
of the water waves problem, we construct an approximate solution as follows : 

• We first construct what we call here approximate initial data, by applying the two successive changes 
of variable on the data (Vq , t]q) : 

f ^,0 - (1 + W - b )) I 1 + K 1 ~ ° 2 ) A ) Vo , 
1 ^s.o = Vo ■ 

• We then choose the parameters (0, A, fi) € [0,1] x R 2 such that the system (T e £ ) belongs to the 
class S of completely symmetric systems (this choice is always possible as we saw previously). Using 
Proposition ^. 6i we know that there exists a unique solution to this system with initial data (V£ , 77^ ) 
: we denote this solution by (V^,rj^). 

• From this exact solution of the symmetric system (T s ^„), we finally construct an approximate solution 
of the water waves problem by successively and approximatively inverting the two changes of variable 
as shown below : 



app 

'fapp = ? /s 



This formal construction of an approximate solution founds its mathematical justification in the following 
theorem which is the last result of this section. 
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Theorem 2.7. Let Ti > 0, s > f + 1, a > s + 3 and (V$j,»7o) be in H a {W l ) d+1 . Let (ip E , rf) Q<£<eo be a 
family of solutions of |0.7p with initial data (V^gi^o) and such that (V?/> £ , 77 e )o<e<e * s bounded with respect 
to e in W 1 - 00 ^, 21]; H a (R d ) d+1 ). We define V s = Vt[ e and choose (9,X,fi) G [0, 1] x M 2 such that the system 
(Tg G E. Then for all e < Eq, there exists T <T\ such that we have : 

Vi £ [0, j] , \V S - V^ p \ L ^ [m , Hs) + |rf - i^pU-do,*]^.) < CeH . 



Proof. We follow in this proof the strategy put forward in [6] : estimates are done on the symmetric system that 
provides the approximate solution rather than on the initial system (|0 . T[) . 



To this end, we take (9,X,fi) G [0, 1] x W 2 such that the system Tg ^ „ belongs to the class £ of completely 
symmetric systems. 

Since (ip s , r] 6 )o< £ <e is a family of solutions of (|0.7j) such that (V?/; e , f? e )o<£<e i s bounded with respect to e in 
W 1,oc ([0, 21];iP(]R d ) d+1 ), using Proposition O implies that (V/, 77 e )o<e< £o wher e = (l + §(1 - 2 )A) V e 
is consistent with the system S g \ . 

Moreover, Proposition 12.41 states that any family (V s ,r) £ )o <E<eo consistent with the system S g I is, up to 
the aforementionned nonlinear change of variables, consistent with the system T e jj . Applying this result to 
(Vg, r? E )o<e<e shows that the family (Vg , V £ )o<s<e i where Vg = (l + |(?7 e — 6)) V e £ , is actually consistent 
with the symmetric system T g I 

Thanks to Proposition ^. 61 we know that there exists a time To such that there exists a unique solution (V£, 77^) 
to this system with initial data (V'so'^so) (defined in the previous formal construction of the approximate 

solution). We are now interested in computing the error estimates between (Vg E ,rj s ) and (V£,r)§,). To this end 
we define V = Vg — V£ and r\ = rf — Tj§^. Writing the equations satisfied by V and 77 and performing standard 
energy estimates on it leads to the following estimate : 

Vi G [0, -±] , |f |£ac ([0 , t] .tf.) + |j?|L-([o,t];H.) < CeH 

where T = min(To,Ti). Inverting the nonlinear change of variables and the pseudo-differential one yields the 
final result. □ 



Remarks 2.8. 

• The construction of the approximated solution of the water waves problem relies on the choice of the 
three parameters 0, A, ^ such that the system Tg\ is completely symmetric. A great advantage of this 
method is that this choice is totally free : we are indeed allowed to choose any suitable triplet {6, A, /1) 
we want, and contruct our approximate solution from the exact solution of the system Tg ^ . In other 
words, approximate solutions of the water waves problem can be constructed starting from the exact 
solution of any symmetric system of the class S. 

• Our theorem relies implicitly on the existence of a family (ip E , V £ )o<e<e °f solutions to the water 
waves problem in Sobolev spaces : in ID-surface, this existence have been already proved by Craig [9], 
Schneider- Wayne [24] thus this implicit hypothesis of existence of solutions is actually a fact. However, 
in 2D-surface, we have no existence result for the water waves problem on a long time scale. Lannes 
proved reccenlty in [16] the existence of solutions to this problem in Sobolev spaces in 2D-surface, but 
we do not know if these solutions persist on a long time scale. Consequently, the analysis is not totally 
complete in 2D-surface. 
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3. The regime of strong topography variations 



In this section, attention is given to the regime of strong variations of the bottom topography : these variations 
are here of order O(l). We recall the Boussinesq-like system (23 2 ) derived in the previous section, and the fact 
that solutions of the water waves problem are consistent with this system. 



Like in the previous section, we aim at deriving asymptotic models, constructing approximate solutions of the 
water waves problem, and justifying these approximations. However, the method introduced in the previous 
section cannot be applyied in the exact same way : this regime is indeed much more complex since the bottom 
terms have here a greater influence than in the first regime. These bottom terms introduce new difficulties 
which compell us to revise and adapt our strategy. 

3.1. A first equivalent system 

First remark that the bottom term h (recall that h = 1 — b is the non-dimensional still water depth) appears 
in the first order term of the second equation of (B2) whereas it is not present in the first one : we have here 
a dissymetry of the order one terms. This fact becomes important when it comes to the BBM trick which is 
unlikely to symetrize these terms. To correctly deal with this regime, we have to invert the order of the change 
of variables, and proceed with an adapted nonlinear change of variables first that symetrizes both order one 
terms and non-linear terms. 

Taking into account the fact that we have to symmetrize both order one terms and nonlinear terms, we introduce 
the following change of variables : 



Assuming that V x V = 0(e), we formally derive the following system of equations satisfied by V and i] : 



(B2) 



d t V + \7 v +^\7\V\ 2 = , 

dtv + V • (h V) + e [V • (r]V) - ±V • VV • V - h 2 \7\7 ■ (h K)) 



= 




so that 




< 



' d t V + VhVr, + ^= [±Vr? 2 + ±V\V\ 2 + (V ■ V)V + VV ■ V 
+ i (i(Vh-V)V-\V\*Vh)] =0(e 2 ) , 

d tV + v(Vh ■ v) + ^ [v • ( V v) - VKv ■ (^vv • (^) 




We introduce the system (7^) that corresponds to the homogeneous version of the previous system 
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where 



d t V + Vh\7r) + %F h 
d t 7 1 + V(Vh-V) + § 



V 
V 



V 
V 



V • (^VV ■ - h 2 VV ■ (VhV)^ 



F,, 

fh 



V 

V 
V 



7^ (1V77 2 + ±V|^| 2 + (V ■ V)V + VV ■ V + i (i(V/i • V)V - \V\ 2 Vh) ) 

A(v-0/V)-&vfc.v) • 



On this new system (7b), we have the following result of consistency : 

Proposition 3.1. Consider a family (ip e , ?y e )o<£<e of solutions of l{0.7\ ) such that (Vip E , T] 6 )o <e<£o is bounded 
with respect to e in W /1:OO ([0, -j-]; H a ($L d ) d+1 ) with a large enough. Then the family {V e ,rf)a <e<eo is consistent 
with the system (7~b), where V e = Vi/ ,£ - 



Proof. First remark that since the velocity held V s is irrotationnal, we have V x V e — 0(e). And since 
(Vt/) 6 , ?7 e )o<£<e i s consistent with the Boussinesq-like system (B2), the previous computations yield directly the 
result. □ 



3.2. Derivation of a class of equivalent systems 

In the previous section, we saw that a suitable change of variable comes from considering Ve, the horizontal 
component of the velocity at the height —1 + 8 (6 £ [0, 1]), instead of the horizontal component of the velocity 
held at the free surface. We can remark that the link between these two variables (and hence the adequate 
change of variables) can be derived from the expression of u app computed during the asymptotic expansion 
process of the operator Z e (er], [3b), which implies that we must adapt our change of variable for our regime of 
strong variations since the expression of u app relies on the considered regime of bottom topography. 

Indeed, we saw in the previous section that the computation of the asymptotic developpment of Z e (j], b)ip relies 
on finding an approximate solution of the elliptic problem (H) on the band 5= [-l,0]xK 2 . Starting from the 
truncation of the computed value of u app at the order 0(s 2 ), 



{z + lf 



h 2 Aip - zhVh ■ Vij) 



+ 0(e 2 ) 



where t/> is the value of the velocity potential at the free surface, shows that Vu app (-, z) gives an approximation 
at order e 2 of the horizontal component of the velocity held, namely V(-, z) — V (/>(■, z) at height z £ [—1, 0]. 
Consequently, in presence of huge bottom variations, the adequate change of variables is given by : 



Vn 



i-U - i)(0v(h 2 v- ) + vv-(/i 2 . )) 



V 



so that 

V = [l + |(0 - l)(9V(h 2 V- ) + VV-{h 2 . ))] V e + 0(e 2 ) 



26 



TITLE WILL BE SET BY THE PUBLISHER. 



From this change of variables, we easily compute the expressions of dtV and V • VhV which we plug into the 
system (%). By rewriting carefully the bottom terms in order to make the quantity \fh appear, one gets the 
following system : 



d tV + W{Vh-V e ) + l 



fh 



v((0 2 - l)h 2 V ■ d t Vg + 2(9 - l)W/i ■ d t V e )] = 0(e 2 ) , 
- V • (> 2 - i)ft 2 VV • {VhVe) + (| 9 2 - |)WW ■ (VhVg) 



-^^VhVh(Vh ■ V ) - - 26 + |)/iVW(V/i • V&)) ] - 0(e 2 ) . 



At this point a new problem arises. Applying the BBM trick in the exact same way as in the previous regime 
leads to a system that is never symmetric for any values of the parameters 6>, A and /i. Indeed, it implies to 
solve a numerical system on the unknowns 9, A and \i which is over-determined. To deal with this problem, 
we simply introduce an additionnal unknown during the BBM trick process : we remark that the term dtVe 
appears twice in the dispersive terms of the first equation, we can then use two different expressions of dtVe, 
each with a different unknowns. This process is summed up in the following relations where we introduce the 
parameters Ai, A2 and fi : 

{ d t Vg = (1 - XJdtVe - \1VhVr1 + 0(e) , 
d t Vg = (1 - \ 2 )d t V e - \ 2 VhV V + 0(e) , 
, V • (VhVg) = AiV • (VhVg) - (1 - fi)dtv + 0(e) , 

where we use the first relation on the term (0 2 — l)h 2 \7-d t Ve and the second relation on the term 2(6— l)hVh-d t V$. 
The last difficulty stands in the possibility to lose the possibly symmetric structure during the derivation of 
the final system. The key point is to rely on the quantity \fh which is crucial to write a symetric form of the 
dispersive terms. 

Finally, we formally derive a new class (Sb) of systems, and we can prove that if a family (V £ 1 rj e ) a<E<£a is 
consistent with the system r fc then (V e e , ?f )o< e < eo where V g £ = [l - §(0 - l)(W(/i 2 V- ) + VV • (h 2 . ))] V £ is 
consistent with any of the following systems (Se,x 1 ,\ 2 ,^) ■ 



(SeMM-.n) < 



(I - ivl) d t V + VhVv 



{l-iV 2 )d t r 1 + V(Vh-V) + 



V 



fh 



V 



+ hVhNV ■ (h 2 V v ) + b 2 s/hV(hVh ■ Vr,) 
+63V/1V • (hVhVv) + b 4 VhVh(Vh ■ Vrj) 



= , 



+ ciV • (/i 2 VV • (VhV) + c 2 V • (WW • (VhV)) 



+c 3 V • (hVhV(Vh ■ V)) + c 4 V • (VhVh(Vh ■ V))) 

where the operators and V\ are defined by 

' V\ = (1-0)((1- Ai)(^ + l)V(/i 2 V- ) + 2(1 - A 2 )V(W/i- )) , 
V 2 = (l-n)((6 2 -\)V-(h 2 V ) + (l9 2 -l)V-(hVhx )) , 



. 



TITLE WILL BE SET BY THE PUBLISHER 



27 



and the parameters (<ii)i<i<4, (6;)i<i<4 have the following expressions : 

f 6i = Ai(l-# 2 ); ci=^ 2 -|); 

b 2 = (1 - #)(2A 2 - |Ai(l + 0)); c 2 - M|^ 2 - |); 

6 3 = ^(l-0 2 ); c 3 = -|0 2 + 20-|; 

I 6 4 = (l-0)(A 2 -^(l + 0)); c 4 = ±(#-2) 2 ; 

The previous computations are summed up in the following proposition. 

Proposition 3.2. Let 9 € [0,1] and (ip £ ,T] £ )o <E<eo a family of solutions of $0.7\ ) such that (V^ e , r] e ) 0<e<Eo 
is bounded with respect to e in W l '°°(\fl, —]; H a (R d ) d+1 ) with a large enough. We define V e = V-0 £ and 

V e = (1 - §(6> - l)(6»V(/i 2 V- ) + VV • (h 2 . ))) (\A;+§^) U e . Then for all (Ai,A 2 ,/x) € R 3 , the family 

(V e ,n e )o <£<£o is consistent with the system (S , 0.Ai.a 2 ,m)- 

Moreover, we have the following proposition on the existence of a subclass of (Sb) composed with fully symmetric 
systems. 

Proposition 3.3. There exists at least one value of (9, Ai, A 2 , /i) smc/i i/iai i/ie system (5 i 6/ ! Ai,a 2 ,m) * s 
symmetric. 

Proof. We are concerned here with the resolution of the following system : 



h 


= ci , 


l>2 


= -c 2 


h 


= c 3 , 


b± 


= ~Ci 



This system on (9, Ai, A 2 , fi) have at least one solution that gives the following approximate values : 

9 w 0.6318 , 
Ai « -0.3416 , 
A 2 w -2.8209 , 
/xw -3.1157, 

which ends the proof. □ 
From now on, we only consider this solution and its approximate values. 

3.3. The fully symmetric systems 

Thanks to Proposition 13.31 we know that some of the systems (Se.\ lt \ 2 .^) of the class (A;,) are completely 
symmetric : we hence denote by £& the non-empty subclass of (At) composed with these symmetric systems. 
Unfortunately, we do not have the same kind of existence theory on these systems as in the previous regime. 

Indeed, the main difference consists in the order one terms of the two equations ( ^f^-p V In order to 



focus on the problem, we rewrite these terms : A(X,dx) I J where A(X, dx ) = ( y ^ 



V • (VhV) 

V \ , A , v Q , f VhVx 
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The proof of the existence of solutions on a short time scale is not modified by these terms, the classical proof 
is still valid. However, the fact that the matrix A depends on the bottom term h is a real problem as far as the 
long time existence is concerned : indeed, one crucial point of the proof here relies on the size of the quantity 
^ on which we have no piece of information. The only case wherein we are surely able to demonstrate the long 
time existence is the case where V/i is of order 0(e) : the term ^ is then of order 0(1) and we can conclude. 
In all other cases, the classical proof fails to provide a rigourous demonstration of the long time existence of 
solutions to these symmetric systems. Nevertheless, we are able to enounce the following proposition : 

Proposition 3.4. Let s > ^ + 1 and (6, X\, X2, fj) be such that the system (5 i 9 j a 1 ,a 2 ./j) belongs to the class E^. 
Then for all (Vo,f]o) € H s (M. d ) d+1 , there exists a time To independant of e and a unique solution (V,n) E 
C{[0,To};H s (R d ) d+1 )n ^([O, T ]; H s ~ 3 (R d ) d+1 ) to the system (Se,x u x aill ) such that (V,rj)\ t=0 = (V ,r) ). 
Furthermore, this unique solution is bounded independently of e in the following sense : there exists a constant 
Co independent of e such that for all k verifying s — 3k > | + 1, we have : 

\(V,T])\ W k,a O ^ 0Ta yj 1 s-3k ( J Si d ) d+l) < Cq . 

Besides, if we suppose that V/i = 0(e), the previous result becomes valid on the long time interval [0, 



( Vl \ 

Proof. The key point of the proof is to demonstrate that the elliptic operator 1 — | I ^2 J is a positive one. 



We first focus on Vl 



(1 - E -VlV, V) = \V\ 2 2 + |(1 - 9 2 )(l - Ai)|W • V\ 2 2 + e(l - 6)(l - X 2 )(Vh ■ V, h\7 ■ V) 
Using the following inequality (satisfied for all a G R) : 



a 2 „ 1 



(Vh-V,hV-V) <^\hV -V\l + ^\Vh-V\l , 



and taking a 2 = ^ 1+ ^J^ 2 leads to 



Vh-V\ 2 



Using the classical Cauchy- Schwartz inequality leads finally to : 



(1 + 0)(1-Ai) 

2(1+9)(1 
(1-0)(1-A 2 )*|V^ 



2 

2 i 



At this point, if we take a small enough value of e, f.e. e < (T^y^i^Plvfcp ' ^ ensures the global positivity of 
V\. On V\, we use the same method : 



(1 - e -V 2 h n,ri) = HI + |(1 - v)(8 2 \)\hVn\l + |(1 - M )(^ 2 - 7 -)(nVh, hV v ) 
Using the same ideas as previously, one gets : 



(i->,,)s(i-| a^zig | V >|j) M 



2 

2 5 
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Once more, if we take f.e. e < ( ua fl2 _^ 2 ._ , , 2 , we have the global positivity of Vf . 

Consequently, taking e < min( (1 J^t-l^v^ ■ (i- M )(|8 2 -|) 2 |Vfe|| ) ensures that tne operator 1- § f ^| ^ is 
positive. 

At this point, using this result and performing usual energy estimates on the system proves the existence of a 
time T such that there exists an unique solution (V, rj) G C([0, T}; H s (R d ) d+1 ) n C x ([0, T]; iP _3 (R d ) d+1 ) to the 
system. □ 

This result gives us an efficient theoretical background to contruct approximate solutions of the water waves 
problem on a time scale 0(1), and O(-) in the case Vh = 0(e). 

This contruction follows the same steps - but in a different order - as the contruction of approximate solutions 
for the first regime : we consider a solution (ip £ ,rj £ ) to the formulation (|0.7| of the water waves problem. We 
take initial data (V'oi^o) sucn that (^V'o'^o) e H s (R d ) d+1 for a suitably large value of s. We then define 
V 6 = V-0 £ and Vq = V^q : we first contruct the data (V£ ,r]£ ) by applying the two successive changes of 
variable on the data (Vq,t]q). We then choose the parameters (9, Ai, A2/i) G [0, 1] x R 2 such that the system 
(Tg ^ Aa ) is completely symmetric. Using Proposition 13.41 we know that there exists a unique solution to 
this system with initial data (V£ ) : we denote this solution by (V£, ryf,). From this exact solution of 
the symmetric system (Sg t \ lt x 2 ,n)t we finally construct an approximate solution of the water waves problem by 
successively and approximatively inverting the two changes of variable as shown below (which is possible if e is 
small enough) : 

VUv = (^-l^^ + l^-^^V^V-^+VV-^FI))) 

Vapp = vh 
We are now able to enounce our final result : 

Theorem 3.5. Let T x > 0, s > f + 1, a > s + 3 and (VV>o,»7o) be in H°(R d ) d+1 . Let (*l> £ ,n £ ) 0<e<£o be a 
family of solutions of |0.7| j with initial data (V^qj^o) an d su ch that (\7ip £ , r/ £ )o <£<eo * s bounded with respect to 
e in W 1 ' oc ([0,T 1 };H a {R J ) d+1 ). We define V £ = Vip £ and choose (0,Ai,A 2 /i) G [0, 1] x R 2 such that the system 

Then for all e < Eq, there exists a time T <T% such that we have : 

\V £ - V app\L^([0,T];H s ) + I 7 ?" - Vlpp\L°°([0,T\;H') < C E 2 

Besides, if we suppose that \7h = 0(e) then {V £ pp ,r) £ pp ) approximates the water waves solutions on a long time 
scale : 

vt g [o, -] , \v £ - ^U-ao,*]^-) + W - ifa P p\L~({o,t];Hn < CeH 



Proof. The proof is an adaptation of the one of Theorem 12.71 and we omit it here. □ 



Remark 3.6. In the general case, where we have no piece of information on the size of the quantity — , our 
analysis is complete on a short time scale. We have indeed an approximation on this interval of time, and we 
know from Lannes [16] the existence of solutions to the water waves problem on a short time scale in 2D and 
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3D. However, in the case V/i = 0(e), this analysis is only complete in 2D - like in the first regime - since we do 
not know about the existence of solutions to the water waves problem on a long time scale. 



This work was supported by the ACT Jeunes chercheurs du ministere de la Recherche "Dispersion et nonlinearites" . 
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